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We discuss an effective theory for QCD at finite chemical potential 
and non-zero temperature, where QCD is reduced to its center degrees 
of freedom. The effective action can be mapped to a flux representation, 
where the complex phase problem is solved and the theory accessible to 
Monte Carlo techniques. In this work, we use a generalized Prokof'ev- 
Svistunov worm algorithm to perform the simulations and determine the 
phase diagram as a function of temperature, quark mass and chemical 
potential. It turns out that the transition is qualitatively as expected for 
QCD. 
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1. Introduction 

Currently one of the main goals of particle physics is to gain a deeper 
understanding of the QCD phase diagram. In the next years running and 
forthcoming heavy-ion colliders will provide access to the regime of QCD 
phase transitions and challenge theory to describe qualitatively and quan- 
titatively the phase diagram structure. Due to the non-perturbative nature 
of QCD, particularly in the transition region, Lattice QCD turns out to 
be the most suitable tool to describe it. However, at finite density Lattice 
QCD faces with the complex phase problem, where the Boltzmann factor 
becomes complex, making the theory inaccessible to Monte Carlo simula- 
tions. This is the main reason for the very limited progress of Lattice QCD 
on the description of the QCD phase diagram. So far the region of small 
chemical potential is accessible to Lattice QCD [Ij, while the rest of the 
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phase diagram is still unexplored. Therefore, it is necessary to develop new 
ideas and methods, which as a preliminary step can be tested in effective 
theories. 

For quenched QCD the deconfinement transition is related to the center 
group ^3 of SU{3) [2j. When the dynamics of the quark fields is coupled, 
center symmetry is broken explicitly, and for /x = it is well established 
that the transition is a crossover [3]. However, one may expect that the 
underlying symmetry still governs parts of the dynamics of the full theory 
[3]. Therefore, in order to study the role of center symmetry for the QCD 
phase diagram, we analyze an effective theory which contains the leading 
center symmetric and center symmetry breaking terms. This effective center 
theory can be mapped exactly to a flux representation [5j, where the complex 
phase problem is solved, and which can be simulated using a generalized 
Prokof'ev-Svistunov worm algorithm [6]. Here we describe the model, the 
algorithm and the results 

2. The Effective Center Theory 

The effective center theory is defined by the action 

X \ i/ = l 

where the dynamical degrees of freedom are the elements of the center group, 
Px G ^3 = {1, e*^'^/^, e"*^'^/'^} at the sites x of a 3-dimensional lattice. 
Px plays the role of the local Polyakov loop, which in quenched QCD is 
the order parameter for center symmetry and thus for confinement. The 
partition function is a sum over all configurations of the center variables, 

z = E{P}eM-s[P]). 

The nearest neighbor interaction term in ([1]), which is invariant under 
global center transformations, may be obtained from a strong coupling ex- 
pansion. The parameter r is an increasing function of the temperature T of 
the underlying Lattice QCD theory. The second term of the effective action 
([TJ may be obtained from a hopping expansion of the fermion determinant 
and contains the leading center symmetry breaking contributions. The pa- 
rameter fi is the chemical potential, and k is a decreasing function of the 
quark mass m. 

For vanishing k, i.e. quenched QCD, the model reduces to the 3d 
3-state Potts model, which is known to have a first order transition at 
r = 0.183522(3) [8]. Increasing k at ^ = the transition weakens until 
it terminates in a critical end point at (r, k) = (0.183127(7), 0.00026(3)) 
[S]. For small non-zero the chemical potential mildens the transition and 
shifts the critical endpoint towards smaller values of k. [51 110 [ [TT \ I12|. 
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3. Flux representation 

Using high temperature expansion techniques one may map the partition 
function to a flux representation ^5j. Its final form is given by: 



{b,s} \x,u / \ X / X \ u / 

(2) 

Here, C = (e^^ + 2e-^)/3, B = (e^^ - e-^)/3C, and the monomer 
weights Ms for s = —1, 0, +1 are: 



J 2k cosh fi 



+ 26"^^°'^'^'^ COS (^\/3KSinh;U - yS)^ 



(3) 



In the new representation ^ the configurations are closed paths formed 
by fiuxes, b^^u £ { — 1,0, +1} which live on the links (x,!^), and monomers 
Sx € {— 1,0, +1} attached to the sites x. They must obey the constraint 
given by the second factor in ([2|). The triality function T(n) = (5o,n mod 3 
enforces the total sum of fiuxes and monomers to be a multiple of 3 at each 
site X. Each configuration of dimers and monomers comes with a strictly 
positive weight and thus the complex phase problem of this model is solved. 

We focus on bulk observables such as the expectation value of the 
Polyakov loop P, its corresponding susceptibility xPj the internal energy 
U and the heat capacity C, which all are obtained as derivatives of the free 
energy. In the end all our observables are expressed in terms of the total 
number of fiuxes and monomers. 



4. Simulation with the worm algorithm 

The most suitable algorithm to perform the Monte Carlo simulation of 
our model is a generalized form of the Prokof 'ev-Svistunov worm algorithm 
[6] . Each worm configuration is generated using four different steps (Fig. [1]) . 
The worm starts at a random position of the lattice (1). It may decide 
to insert fiuxes and move to the neighboring site (2) but can also insert 
monomers (3). The insertion of a monomer is followed by a random hop to 
another position (3), where again a monomer is inserted. These steps are 
continued until the worm closes (4). 

The results from the new algorithm were checked using several strategies. 
For the quenched case (k = 0) the known results were reproduced For 
small values of r we determined the power series perturbatively for the 
partition sum, taking into account the terms up to r^. We found excellent 
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agreement between the Monte Carlo results and the perturbative series (see, 
Fig. [2]). Finally, two independent programs were written for cross checks. 

We generated ensembles on 36^, 48'^, 64^ and 72'^ lattices, for k = 
0.1, 0.01, 0.005 and 0.001. The evaluation of our observables (P), xp, 
U and C is based on up to 10^ configurations which are separated by 10 
worms for decor relation. Autocorrelation times were determined for several 
observables and used in the estimate for the statistical errors. 



The Ihs. plot of Fig. Oshows (P) /V as a function of r and fi at k = 0.01. 
For small r and ^, {P) /V is rather small and center symmetry is broken 
only very mildly, which implies that matter is confined. On the other hand, 
when r and // increase, the system undergoes a change and {P) /V reaches 
values close to 1. For QCD this means that the system is driven into the 
deconfined phase signaled by a rapid increase of the Polyakov loop. 

To identify the phase boundary in the r-^u plane (Fig. [2]rhs. plot), we 
use the position of the maxima of xp- Foi' this determination we fit the 
data for Xp near the maxima with a parabola and obtained the position of 
the maximum as one of the fit parameters. The corresponding statistical 
error was computed with the jackknife method. The points with horizontal 
(vertical) errors were determined at fix r (^u) as a function of fj, (r). We 
compare the results for 4 values of n. The dashed horizontal line at the top 
marks the value of the critical r for n = 0. The dashed curves near the 
bottom of the plot are the results from the perturbative series for small r, 
which we briefly discussed in the last section. 

To determine the nature of the phase transition we compared the phase 
boundaries of the two order parameters and the heights of their fluctuations 
at different volumes. We found that the transition lines from xp C 
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Fig. 1. Schematic illustration of the worm algorithm. 



5. Results from the Monte Carlo simulation 
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Fig. 2. Left: Order parameter {P)/V as a function of r and ioi V = 36'^ at 
K = 0.01. Right: Phase diagram obtained from the maxima of the Polyakov loop 
susceptibihty for 4 values of k. The dashed horizontal line marks the critical value 
of T for K = 0, and the dashed curves at the bottom are the results from the 
perturbative t expansion. 




Fig. 3. Left: The critical parameters are determined by locating the position of 
the maxima of the susceptibility xp- Right: Comparison of the phase boundaries 
obtained from the maxima of susceptibility xp and heat capacity C. 



do not coincide (Fig. [3] rhs. plot), and there is no volume scaling (Fig. [3] 
Ihs. plot). This is a clear signal that the phase changes are crossovers. 
Although there is no critical end point when only center degrees of freedom 
are considered, the qualitative behavior and the mass dependence of the 
phase boundaries are as one expects for full QCD. 
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6. Conclusions 

We have studied an effective tlieory of QCD witli finite cfiemical poten- 
tial at non zero temperature. Mapping this theory to the flux representation 
enables us not only to have a model free of the complex phase problem but 
also to generate configurations efficiently in a wide range of parameters us- 
ing a worm algorithm. For all parameter values we studied the transition is 
of a crossover type and we conclude that center symmetry alone does not 
provide a mechanism for first order behavior in the QCD phase diagram. 

Currently we consider two possible further directions: The effective the- 
ory can be made more realistic by replacing the Z3 spins by continuous 
SU{3) valued variables |13j . Also, we have implemented a parallel ver- 
sion of the worm algorithm, to take advantage of the compute capability of 
graphic cards. 
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